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Abstract 

We calculate explicitly the Bethe vectors states by the algebraic Bethe ansatz method 
with the 5fZ(2)-invariant i?-matrix for the two-site Bose-Hubbard model. Using a binomial 
expansion of the n-th power of a sum of two operators we get and solve a recursion equation. 
We calculate the scalar product and the norm of the Bethe vectors states. The form factors 
of the imbalance current operator are also computed. 


1 Introduction 


The first experimental verification of the Bose-Einstein condensation (BEC) [U E] |3] occurred 
after a gap of more than seven decades following its theoretical prediction mn- After its 
realization a great deal of progress has taken place both in the theoretical and experimental 
study of this physical phenomenon [6l[7ll8ll9l[ini[IIl[l2l[l3]. A particularly fruitful instru¬ 
ment in relation to ultracold physics are many-body atomic models related to BEC. In this 
direction the quantum inverse scattering method (QISM) [Ml |T5l [161HTJ [TH] been used to 
solve and study some prototypical many-body models that contribute to describe phenom¬ 
ena associated to BEC UniEolEl]. Some of these models, despite their simplicity, display 
a rich structure showing quantum phase transitions and interesting semi-classical behaviour 
that have been studied in [22l [23l EH EH EH Elj, and explored in different areas such as nu¬ 
clear physics, condensed matter and atomic-molecular physics. To keep things as simple as 
possible we shall consider here the two-site Bose-Hubbard, also known in special cases as the 
canonical Josephson Hamiltonian [7] . This model may be viewed as a particular case of the 
bosonic multi-state two-well model studied in [21], and can be used to describe a quadrupolar 
nuclei system in nuclear magnetic resonance [28] by a 77/2 Schwinger pseudo-spin realization 
of the Hamiltonian. Conversely, there is also a link with the Schwinger bosonic realization 
of the Lipkin-Meshkov-Glick model |29l [30] used to study closed shells in a nuclei model. 
In spite, of course, from being a two-site specialization of the Bose-Hubbard model it is a 
very useful model in various realms as the understanding of tunnelling phenomena using two 
BEC [m |32l |33l Ell ESI ESI ET] , as well as quantum phase transitions using tools of quantum 
computation and quantum information. The model is described by the Hamiltonian 

H = —{Ni — N 2 Y - ~ ^ 2 ) —+ h^Ui), (1.1) 

where, a\,dl, denote the single-particle creation operators in each site and, Ni = d\di,N 2 = 
0 , 20 , 2 , are the corresponding boson number operators. The total boson particles number 
operator, N = Ni + N 2 , is a conserved quantity, [H,N] = 0. The coupling K provides the 
interaction strength between the bosons and is proportional to the s-wave scattering length, 
A/i is the external potential and Sj is the amplitude of tunnelling. 

The Hamiltonian fll.ip is integrable in the sense that it can be solved by the quantum 
inverse scattering method (QISM) and it has been discussed in different ways using this 
method [371 EH ESI SSI IHl SSI SSI SH US]- The algebraic formulation of the Bethe ansatz, 
associated to the QISM, was primarily developed in [TH [151 [Ml EZl [T8] . 
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A very important problem in the algebraic Bethe ansatz method is the construction of the 
Bethe vectors states (BVS) [38l 061 07] using the correspondent creation operator applied 
to the pseudo-vacuum. Employing this form of the BVS it is possible m to calculate their 
scalar product and then use it to calculate important physical quantities as the form factors. 
Form factors are dehned as the matrix entries of operators in the base of the eigenvectors of 
the Hamiltonian. Another important application is in the calculation of the average values 
of the operators as for example correlation operators. Applying this method, some physical 
quantities for the Hamiltonian (II.ip were obtained in [3H]. Recently, BVS have shown to be 
useful in fundamental issues of planar A/" = 4 super Yang-Mills (SYM) theory [08109] in the 
context of the integrability in the AdS/CFT correspondence [50]. 

In the present work, we develop and use a new method to explicitly calculate the BVS and 
obtain the scalar product of two BVS for the two-site Bose-Hubbard model. Although we 
concentrate on this model the procedure is of general applicability. We use a Lax operator 
to construct a realization of the monodromy matrix and get an algebraic identity between 
the associated C-operator and the D-operator of the monodromy matrix (see next section) 
needed to calculate the BVS. We then use the binomial expansion for the n-th power of 
the sum of two operators: in a hrst step we will show that the binomial expansion can be 
written as a sum of permutations of the product of that operators or as a standard binomial 
expansion as in a commutative algebra plus a function of the commutator of these two 
operators; in a second step we write a recursion equation and give its solution. Next, we 
calculate the scalar product between one on-shell and one generic off-shell BVS, as well as 
the norm. As an application we obtain the form factors (non normalized) for the imbalance 
current operator. 
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2 The algebraic Bethe ansatz method 


The spectrum of the Hamiltonian (II.ip has appeared in different papers [371 EHl EHl HOI 021 
031 mi 05] using this method. To £x notation we will shortly describe the algebraic Bethe 
ansatz method, see [391001 [51] for more details. We begin with the 5fZ(2)-invariant i?-matrix, 
depending on the spectral parameter u, 


/ 1 0 0 0 \ 

^ 0 b{u) c{u) 0 

0 c{u) b{u) 0 ’ 

\ 0 0 0 1 / 

with b{u) = u/{u + 7]), c{u) = rj/{u + rj) and b{u) + c{u) = 1. 
parameter, to be chosen later. 

It is easy to check that R{u) satisfies the Yang-Baxter equation 


( 2 . 2 ) 


Above, Tj is an arbitrary 


Ri 2 {u - v)Ri3{u)R23{v) = R23{v)Ris^u)R^i^u - v), (2.3) 


where Rjk{u) denotes the matrix acting non-trivially on the j-th and the fc-th spaces and as 
the identity on the remaining space. 

Next we define the monodromy matrix T{u), 


f{u) 


A{u) B{u) \ 

C{u) D{u) ) ’ 


such that the Yang-Baxter algebra is satisfied 


(2.4) 


Ri 2 {u - v)fi{u)f 2 {v) = f 2 {v)fi{u)Ri 2 {u - v). (2.5) 

To obtain a solution for the two-site Bose-Hubbard model (II.Ih we need to choose a real¬ 
ization for the monodromy matrix 7r(T(n)) = L{u). In this construction, the Lax operator 
L{u) has to satisfy the algebra 
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Ri2{u - v)Li{u)L2{v) = L2{v)Li{u)Ri2{u - v), 


( 2 . 6 ) 


where we use the standard notation. 

Then, dehning the transfer matrix, as usual, through 


i{u) = Tr 7i{T{u)) = ni^Aiu) + D{u)), (2.7) 

it follows from fl2.5p that the transfer matrix commutes for different values of the spectral 
parameter. 


We are using the well known [19] Lax operator, solution of the equation fl2.6l) . 


Li{u) 


fui + r]Ni 

I 4 



*= 1 . 2 , 


( 2 . 8 ) 


for the boson operators a|, a,, and W. These operators obey the standard canonical boson 
commutation rules. 


Using the co-multiplication property of the Lax operators (12.Sh we get the following real¬ 


ization for the monodromy matrix, 

n{T{u)) = Li{u + u)L 2 {u — u), (2.9) 

whose entries are, 

7i{A{u)) = — + 

ti‘^NiN2 — tjuj^Ni — N 2 ) T (2.10) 

ti{B{u)) = {u + u+ riNi)a2 + ri~^ai, (2.11) 

TriCiu)) = {u - u + rjN 2 )a\+ (2.12) 
7r(T>(n)) = a|a 2 + (2.13) 


Hereafter we will use the same symbol for the operators and its respective realization, so we 
dehne ti{0{u)) = 0{u) for any operator in the entries of the monodromy matrix fl2.4l) . 

The parameters of the Hamiltonian fll.ip are all real numbers, K, A/r, Sj G M. The 
parameters in the operators (12.10112.11112.12112.131) can be complex numbers, u, r], uj G C, but 
in this case the transfer matrix is not Hermitian. We will only consider the Hermitian case. 
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We can apply the algebraic Bethe ansatz method, using the Fock vacuum as the pseudo¬ 
vacuum |0) = |0)i ® 10 ) 2 , to hud the BAE, 


N 






Vj -Vj-r] 

- Vj +T]' 


i,j = 


(2.14) 


3 Bethe vectors states 


In the algebraic Bethe ansatz method, the BVS are constructed by the application of the 
(7-operator to the pseudo-vacuum |0), 


N 

I*)=n 1°)' (3.15) 

j=i 

where the are solutions of the BAE fl2.14p . 

Using the H-operator fl2.13p we can write the (7-operator fl2.12p as. 


C(v,) = f,a\+-D, 

where 

fj = Vj - w, 

and 

V = 1] all) = rj a\N2 + al 7 ]~^. 


Now we can write the BVS in the product form as 

N 

i=i 

or in the summation form as 


( 3 . 16 ) 

( 3 . 17 ) 

( 3 . 18 ) 


( 3 . 19 ) 


N 

1^) = (3.20) 

n=0 

with the identification 
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(3.21) 


N 


N 


J^O = Y\fv J^N = 




0 


i=i 


ii j2, • • • j'at=1 • • • /. 

(ji>i2>...>iiv) 


JN 


To explicitly write the BVS (13.191) or fl3.20p we need to expand the powers of the V- 
operator fl3.18p . We thus consider the following n-th power binomial expansion for any two 
operators X and Y, with n> 2, proved by induction, 


n 2!(n-j)! 

{x+Yr = 

j=0 0=1 


n-J) 


n ).Y»-^y^+f(|.Y,y]), 

j=0 ^ 


with 


n 


ni 


(3.22) 


(3.23) 


j ) j!(^-i)!’ 

where the Th-operator stand for the set of all permutations of the X and Y operators, 
not respecting the commutation rule between them, and f([X,y]) is a function of that 
commutation rule. Below we show f([X, y]) for two values of n: 


I - For n = 2: 


II - For n = 3: 


f([x,y]) = [y,x]. 


(3.24) 


f([x,y]) = [y,x2] + [y^x] + [xy,x] + [y,xy]. 


(3.25) 


If the operators commute, f([X,y]) = 0, we clearly get the standard commutative bino¬ 
mial formula. In this case the n-th power of the /^-operator fl2.13p is given by 


D" = 


y 

i=o 




(3.26) 


The n-th power of the £)-operator (13.261) is acting in the pseudo-vacuum as 
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w | 0 ) = r]-^^ | 0 ). 


( 3 . 27 ) 


Redefining the P-operator fl3.18l) as 


V = X + Y, (3.28) 

with 


X = r]d\N2, Y = 1] ^dl, (3.29) 

we get the binomial expansion of the "D-operator fl3.18|) for the n-th power, 

n! 

n j!(n-j)! 

= ( 3 . 30 ) 

j = 0 Q;=1 


Applying the n-th power of the "D-operator (I3.3np to the pseudo-vacuum we get 


= |n-j)2, (3.31) 

j=o 

where the coefficients Cnj satishes the recursion equation 


Cn+i,j = V VJ + C'nj-1 + V ^^/nTTXj Cn,j, (3.32) 

with the condition Cjj = 0. 

The solution of the recursion equation fl3.32p is 


^n,j h 






n-j 


[n- i)\ 


1=0 


n 


(3.33) 


Finally, using the binomial expansion of the P-operator fl3.3ip we can write the BVS 
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flCTl as 


IM/) = 


X 


J-Q \N)i 0 10)5 

N n 

E E J 

n=l j=0 

\N -n + j)i 0 \n-j)2. 


{N -n + j)\ 


(3.34) 


The scalar product between one on-shell and one generic off-shell BVS fl3.34p . IT) and 
|T), with the set of solutions of the BAE (I2.14p and the generic set [52], is 


{T|T) = 


X 


N\ ^0* J-Q 

N n 

^ ^ ^ ^ ^r ^n C*r,r—n+j ^n,j 
r,n=l ^=0 

(A^ — n -I- j)! 


From the scalar product fl3.35|] we can write the norm of the BVS fl3.34p . 


(3.35) 


(>i>i>i>) = 

X 

Because the total number of atoms is a conserved quantity all the BVS (13.341) are eigen¬ 
functions of the V-operator with the same eigenvalue V, and so they are degenerate states 
for this operator. As Ni and N 2 are not conserved quantities the BVS are not eigenfunction 
of these operators but we can still write the form factors of these operators using the BVS 
fl3.34p . For instance, the non normalized form factors of the imbalance current between the 
two BEC is written as 


N n 

N\ l-Tol^ ^ ^ Jq JV Cr,r-n+j Cn,j 

r,n=l j=0 

{N — n + jy. 

\/j\ir-n + jy 


(3.36) 
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^ Nt - No 


N 


X 


N n 

r,n=l jf=0 

2{n- i) 


X Cr j-—n+j ^n,j 


N 


{N — n + j)\ 
A/j!(r - n +j)! 


(3.37) 


4 Summary 

We have explicitly written the Bethe vectors states (BVS) by the algebraic Bethe ansatz 
method using the 5f/(2)-invariant i?-matrix and an algebraic relation between the C and the 
D operators. We use a binomial expansion of the n-th power of the sum of two operators 
obtained from that algebraic relation to write a recursion equation and solve it. The binomial 
expansion of the n-th power of the sum of two operators can be written as a commutative 
binomial expansion plus a function of the commutator of the operators. We calculate the 
scalar product and the norm of those BVS. The BVS are degenerate eigenfunctions of the 
total number of particles iV-operator with eigenvalue N. As an example of application of 
the BVS we calculate the form factors for the imbalance current operator. 
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